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Context. The blazar AO 0235+164 was claimed to show a quasi-periodic behaviour in the radio and optical bands in the past, with the main 
outbursts repeating every 5-6 years. However, the predicted 2004 outburst did not occur, and further analysis suggested a longer time scale, 
j"^ , according to which the next event would have occurred in the 2006-2007 observing season. Moreover, an extra emission component contributing 
to the UV and soft X-ray flux was detected, whose nature is not yet clear. An optical outburst was observed in late 2006 - early 2007, which 
triggered a Whole Earth Blazar Telescope (WEBT) campaign as well as target of opportunity (ToO) observations by the Swift satellite. 
Aims. In this paper, we present the radio-to-optical data taken by the WEBT together with the UV data acquired by the UltraViolet and Optical 
Telescope (UVOT) instrument onboard Swift to investigate both the outburst behaviour at different wavelengths and the nature of the extra emission 
r- H ■ component. 

^ ' Methods. Multifrequency light curves have been assembled with data from 27 observatories; optical and UV fluxes have been cleaned from the 
' contamination of the southern active galactic nucleus (AGN). We have analysed spectral energy distributions at different epochs, corresponding to 
different brightness states; extra absorption by the foreground galaxy has been taken into account. 
£S) ■ Results. We found the optical outburst to be as strong as the big outbursts of the past: starting from late September 2006, a brightness increase of 
7-H ' ~ 5 mag led to the outburst peak in February 19-21, 2007. We also observed an outburst at mm and then at cm wavelengths, with an increasing 
time delay going toward lower frequencies during the rising phase. Cross-correlation analysis indicates that the 1 mm and 37 GHz flux variations 
lagged behind the R-bmd ones by about 3 weeks and 2 months, respectively. These short time delays suggest that the corresponding jet emitting 
regions are only slightly separated and/or misaligned. In contrast, during the outburst decreasing phase the flux faded contemporaneously at all cm 
wavelengths. This abrupt change in the emission behaviour may suggest the presence of some "shutdown" mechanism of intrinsic or geometric 
nature. The behaviour of the UV flux closely follows the optical and near-IR one. By separating the synchrotron and extra component contributions 
to the UV flux, we found that they correlate, which suggests that the two emissions have a common origin. 

Key words, galaxies: active - galaxies: BL Lacertae objects: general - galaxies: BL Lacertae objects: individual: AO 0235+164 - galaxies: jets - 
galaxies: quasars: general 
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1. Introduction data acquired by the XMM-Newton satellite during 3 point- 

ings in January and August 2004, and January 2005, and by 
AO 0235+164 at redshift z = 0.94 is one of the best-studied optical spectroscopic observations with the 3.6 m Telescopio 
BL Lac objects. The analysis of its radio and optical light Nazionale Galileo (TNG ). The results of this inten s e obser ving 



curves extending over 25 years led | Raiteri et al] (EooTb to sug- effort were published bv dRaiteri et al.l20 05. 2006^ 120073: see 



gest a quasi-periodic occurrence of the main outbursts ev- also Hage i^Thorn et al.ll2007bl for an analysis of colour variabil 

ery 5.7 ± 0.5 years. The next outburst was predicted around ity) The pred i c ted outburst was not observed, and time-series 

February-March 2004, and a multiwavelength campaign was analyfjis on the Hght curyes exte nded to 2005 revealed a possible 

organised by the Whole Earth Blazar Telescope (WEBTJJ to c hai-acteristic variability time scale of ~ 8 years. Moreover, the 

follow the expected event. The radio-to-optical observations by XMM-Newton observations suggested the presence of an extra 

the WEBT were complemented by the optical-UV and X-ray emission component in the source spectral energy distribution 

(SED), in addition to the synchrotron and inverse-Compton ones. 

Send offprint requests to: C. M. Raiteri The origin of this component, peaking in the UV/soft X-ray fre- 

* The radio-to-optical data presented in this paper are stored in the quency range, was ascribed either to thermal emission from an 

WEBT archive; for questions regarding their availability, please contact accretion disc, or to synchrotron emission from an inner jet re- 

the WEBT President Massimo Villata (villata@oato.inaf.it). ^ Qn 

1 htt p : //www . to . astro . it/blazars/webt/ 

(see e. g. TviTlata et al., 2006, 2007; Boettcher et al., 20071; iRaiteri et all An increased r adio activity was d etected in the 2005-2006 



:tivity was a e 

2007b) observing season (iBach et all 120071) . and a dramatic optical 
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brightening was finally observed in late 2006 - early 2007. This 
triggered a new WEBT multiwavelength campaign as well as 
ToO observations by the Swift satellite. In this paper, we present 
the radio-to-optical observations performe d from spring 2005 
(i.e. the end of the period considered in iRaiteri et all 12006a) 
to October 2007, and the data acquired by the Ultra Violet and 
Optical Telescope (UVOT) instrument onboard Swift. X-ray 
data from the Swift X-ray Telescope (XRT) and Burst Alert 
Telescope (BAT) instruments will be presented in another paper 
(Kadler et al., in preparation). 

This paper is organised as follows: the procedures we 
adopted to treat the WEBT and UVOT data are outlined in 
Sect. 2, with particular attention paid to the subtraction of the 
southern AGN contribution from the optical and UV photome- 
try. The multifrequency light curves are presented and discussed 
in Sect. 3 while in Sect. 4 time lags among variations at dif- 
ferent frequencies are derived by means of statistical analysis. 
Spectral energy distributions with simultaneous data from near- 
IR to UV are constructed in Sect. 5, with the aim of separating 
the synchrotron and the extra emission components and of un- 
derstanding their relationship. Finally, the results of our work 
are discussed in Sect. 6. 



2. Data reduction and analysis 

2.1. Radio-to-optical observations 

The new WEBT campaign on AO 0235+164 saw a wide inter- 
national participation. Dat a were collec t ed start ing from the end 
of the period considered in lRaiteri et al.l(l2006al) . i.e. from spring 
2005 to October 2007. Table[T]shows the list of participating ob- 
servatories, the size of the telescopes, and the bands/frequencies 
at which we acquired the data. 

We collected optical and near-IR data as instrumental mag- 
nitudes of the source and comparison stars in the same field. 
Magnitude c alibration w a s perf ormed using preferably Stars 1, 
2, and 3 by ISmith et alj dl985l) in the optical bandfl and the 
iGonzalez-Perez et alj d2001l) photometric sequence in the near- 
IR. Optic al light curves we re carefully assembled and cleaned 
(see e.g. Villat a etalll2002l) . 

We collected radio data as calibrated flux densities; no sys- 
tematic offset among datasets provided by different observers 
was evident; data binning was needed in a few cases to re- 
duce the noise. We also included data from the VLA/VLBA 
Polarization Calibration Databas^j 1 (PCD) in our light curves. 

2.2. UV observations by Swift 

The Swift satellite observed the source from January 28 to 
February 26, 2007. The UVOT instrument (iRominget all feOOS) 
acquired data in the U, UWW1, UVM2, and UVW2 niters, which 
were processed with the uvotmaghi st task of the HEASOFT 
6.3 package. Following the recommendations contained in the 
release notes, we used an aperture radius of 5 arcsec in agree- 
ment with the standard photometric aperture defined in the cal- 
ibration files (CALDB updated as July 2007). Background was 
extracted in an annulus centred on the source, with inner and 
outer radii of 8 and 1 8 arcsec, respectively. 



2 In some cases, we also used Stars 9, 10, and 11 by 
Gonzalez-Perez et al. (2001). 

3 http : //www . via . nrao . edu/astro/calib/polar/ 



Table 1. Ground-based observatories participating in this work. 



Observatory 


Tel. size 


Bands 


Radio 




Crimean (RT-22), Ukraine 


22 m 


22, 37 GHz 


Mauna Kea (SMA), USA 


8x6m fl 


0.85, 1 mm 


Medicina, Italy 


32 m 


5, 8, 22 GHz 


Metsahovi, Finland 


14 m 


37 GHz 


Noto, Italy 


32 m 


43 GHz 


UMRAO, USA 


26 m 


5, 8, 14.5 GHz 


Near- 


infrared 




Campo Imperatore, Italy 


110 cm 


J,H,K 


Roque (NOT), Spain 


256 cm 


J,H,K 


Optical 




Abastumani, Georgia 


70 cm 


R 


Crimean, Ukraine 


70 cm 


B,V,R,I 


Greve, Italy 


32 cm 


R 


Hankasalmi, Finland 


40 cm 


V,R,I 


Jakokoski, Finland 


50 cm 


R 


Kitt Peak (MDM), USA 


130 cm 


B,V,R,I 


L'Ampolla, Spain 


36 cm 


R 


Michael Adrian, Germany 


120 cm 


R 


New Millennium, Italy 


36 cm 


R 


Roque (KVA), Spain 


35 cm 


R 


Roque (NOT), Spain 


256 cm 


U,B,V,R,I 


Rozhen, Bulgaria 


50/70 cm 


B,V,R,I 


Sabadell, Spain 


50 cm 


R 


Sobaeksan, South Korea 


61 cm 


V,R,I 


St. Petersburg, Russia 


40 cm 


V,R,I 


Torino, Italy 


105 cm 


R,I 


Toyonaka, Japan 


35 cm 


B,V,R,I 


Tuorla, Finland 


103 cm 


R 


Valle d'Aosta, Italy 


81 cm 


V,R,I 


Yunnan, China 


102 cm 


R 



" Radio interferometer including 8 dishes of 6 m size. 



2.3. Correction for ELISA contribution 

A major point in the analysis of the optical and UV radiation 
coming from AO 0235+164 is to elimina te the contaminatio n 
from the southern AGN, named ELISA bv lRaiteriet all d2005l) . 
Since it lies only 2 arcse c south of the source and it is very 
faint ( Nil sson et all 1 1996b . this object is not resolved from AO 
0235+164 in most optical frames as well as in the UVOT ones. 
Consequently, it gives a contribution to the AO 0235+164 pho- 
tometry, which affects more the blue than the red part of the 
optical-UV spectrum and is stronger wh e n the b lazar is fainter. 
To correct for this effect, IRaiteri et al.l (2005) estimated the 
UBVRI magnitudes of ELISA from images taken at medium- 
size telescopes participating in the 2 003-2004 WEBT campaign. 
More recently. [Raiteri et all d2007al) analysed an ELISA opti- 
cal spectrum taken at the 8 m Very Large Telescope (VLT). 
The UBVRI magnitudes derived from this spectrum confirm the 
previous photometry results and lower the uncertainties to 0.05 
mag; only in the case of the V filter a slightly fainter value ap- 
pears more appropriate. Hence, to subtract the ELISA contribu- 
tion from the optical fluxes, we adopted the following ELISA 
magnitudes: U = 20.80, B = 21 AO, V = 20.95, R = 20.50, and 
/ = 19 .90. Zero-magnitude fluxes were taken from lBessell et al.l 
d!998l) . 

We then estimated the ELISA contribution in the UVW1, 
UVM2, and UW2 filters from th e HST /FOS composite 
spectrum published by Burbidge et al.l dl996Sn . We subtracted 

4 We note that the flux in the UVIV1 band is affected by the contri- 
bution of the C m] emission line, and that the HST/FOS UV spectrum 
nicely reconnects in flux with the VLT optical spectrum. 
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Fig. 1. Optical and near-IR light curves of AO 
0235+164 in the 2006-2007 observing season. 
The optical magnitudes are corrected for the 
ELISA contribution. Maximum and minimum 
brightness levels are highlighted with horizon- 
tal dotted lines. 



these ELISA flux densities from the total flux densities ob- 
tained with uvotmaghist and converted the results into mag- 
nitudes with zero-magnitude fluxes of 3.963, 4.544, 5.400 x 
10" 9 ergcnr 2 s -1 A -1 . UV magnitudes of ELISA are: 20.45, 
20.50, and 20.65 in the UW1, UVM2, and UVW2 bands, re- 
spectively, again with a 0.05 mag uncertainty. 

As in the previous works, the ELISA contamination in the 
near-IR is assumed to be negligible, and null in the mm and cm 
bands. 



3. Light curves behaviour 

The final light curves are shown in Figs.QJ[3] The optical and UV 
magnitudes have been corrected for the ELISA contribution, as 
explained in the previous section. 

Figure Q] displays the optical and near-IR light curves of the 
2006-2007 observing season. Notwithstanding the differences in 
time sampling, the general trend is similar at all wavelengths, as 
expected. At the beginning of the observing season the source 
was in a faint state, but by in late September the source bright- 
ness started to increase. We observed a fast flare to peak at JD 



= 2454077.4 (December 7, 2006). It was characterized by a 2.0 
mag brightening in 4 days in the R band, and by noticea ble in- 
tranight variability (see also Ha gen-Thorn et all l2007al) . After 
the drop of 1.5 mag in 3 days, the source brightness continued 
to rise smoothly, with minor flares superposed. The outburst cul- 
minated on JD = 2454151-3 (February 19-21, 2007). This peak 
is well sampled in the optical bands, but was missed by the near- 
IR observations, which registered the brightest level some days 
later, when the optical flux was already in a decreasing phase. 
We notice that in this observing season the variability amplitude 
in the best-sampled R band was ~ 5 mag. 

Observations by the Swift satellite were triggered in late 
January during the culminating phase of the outburst. The U and 
UV light curves from UVOT are shown in Fig.|2]as blue circles, 
and compared to the 7?-band light curve. In the t/-band panel, 
ground-based data from the Nordic Optical Telescope (NOT) 
were plotted as red diamonds. UVOT data have been binned 
daily, and a lower limit to the error of 0. 1 mag has been adopted 
to take into account statistical and systematic errors. The grey 
line in each panel represents a cubic spline interpolation through 
the 1-day binned /?-band light curve (rescaled according to aver- 
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Fig. 2. UVOT light curves of AO 0235+164 compared to the 
ground-based U data from the NOT (red diamonds) and to the 
/?-band data (bottom panel). All data have been corrected for the 
ELISA contribution. The cubic spline interpolation through the 
1-day binned R -band light curve is shown in each panel, properly 
rescaled (see text). 



age colour indices: U-R = 2.1; Wl-R = 2A;M2-R = 2.3;and 
W2 -R = 2.7). It helps to recognize that the trend of the /?-band 
light curve is well followed by all the UVOT light curves. Also, 
the variability amplitude appears nearly the same in the optical 
and UVOT bands. 

The increase of the optical and near-IR fluxes was followed 
soon after by an increase of the millimetric and then of the centi- 
metric fluxes, with a progressive time delay going toward longer 
wavelengths. Figure [3] shows the radio light cu rves starting 
from th e end of the observing period presented in iRaiteri et al.l 
(2006a), compared to the R-band one (top panel). Optical obser- 
vations in the 2005-2006 season were rather sparse; they show 
moderate source ac tivity at a level sim ilar to that observed in 
the previous season (Raiteri et al, 2006a). The figure also shows 
that, after the solar conjunction in 2007, the source was again in 
a very faint optical state. 

The second panel from the top displays 1 mm and 850 fim 
data from the SubMillimeter Array (SMA). A cubic spline inter- 
polation through the 10-day binned 1 mm light curve is plotted 
to better follow the behaviour of the source flux at this wave- 
length. The 1 mm flux density increased by a factor ~ 8 from 
the beginning of the 2006-2007 season to the time of the ob- 
served maximum, on March 6, 2007. After this date, the only 
data point acquired before the solar conjunction epoch suggests 
that the outburst peak was already over. Moreover, data after so- 



lar conjunction show that the mm flux density had reduced by 
more than a factor of 2. We notice that a first hump in the 1 mm 
light curve, around JD = 2454100, may be correlated to the fast 
optical flare observed at JD ~ 2454080. 

In the following panels, radio light curves at lower and lower 
frequencies are displayed. We have plotted the cubic spline in- 
terpolation through the 20-day binned light curve at 37 GHz, 
the best-sampled one, to allow an easier comparison of the flux 
density behaviour at different wavelengths. In general, the ra- 
dio spectrum appears "inverted": the flux density is higher at 
higher frequencies. This is not true in the period (JD ~ 2453900- 
2454050) leading to the minimum flux levels, where the flux 
density appears higher at longer wavelengths, and during the 
outburst decreasing phase, where the radio spectrum seems to 
be flat. 

The radio emission had already undergone a rather bright 
phase in 2005-2006, a kind of "minor outburst" like those of- 
ten observed in the past, which sometimes are accompanied by 
noticeable optical events and sometimes not. A flux minimum 
was reached around JD = 2454000-2454100, and then the flux 
started to increase again, all these variations occurring first at 
the higher and then at the lower frequencies. However, the 22 
GHz flux density reached a maximum value similar to that ob- 
served at 37 GHz almost at the same time, and the flux decrease 
was also similar and contemporaneous. This is highlighted in the 
fourth panel of Fig. |3]by the overlapping of the spline through 
the 37 GHz data with that through the 22 GHz ones. The lower- 
frequency light curves are more sparsely sampled, but they too 
suggest a progressively delayed rising phase, and a contempora- 
neous decreasing phase. 

While the general multifrequency behaviour during the out- 
burst seems in a greement with the predictions of shock-i n-jet 

aet al.1 1 1992b . the 



models (see e.g. (Hughes et al.L 1 19891 



preaictio 
<■ Valtaoj 



decay stage might suggest a more complex picture. Indeed, it 
seems as if a kind of "shutdown" occurred when the perturbation 
that was responsible for the flux increase reached the jet region 
emitting at cm wavelengths. The shape of the 37 GHz outburst, 
sharp and without a plateau, may indicate that the outburst ob- 
served at this frequency has already been affected by the shut- 
down. The signature of the shutdown would be mainly visible in 
the 22 GHz light curve, which is characterized by an abrupt flux 
decay. 

It is interesting to compare this recent behaviour of the 
source with its past history. Figure [4] displays the light curves in 
the optical R band and at two radio frequencies in the last ~ 12.5 
years, including the previous radio-optical outburst observed in 
1997-1998. We notice that the last outburst had a shorter du- 
ration with respect to the outburst in 1997-1998. Moreover, if 
we compare the 2005-2006 and 2006-2007 events, we see that 
the latter was much "harder" (i.e. stronger at higher frequencies) 
than the former. Similarly, the 1999-2000 post-outburst event 
is harder than the preceding major one. Other very hard (opti- 
cal only) events are visible along all the period shown in Fig. [4] 
starting from the 1995-1996 eve nt dTakalo et all 119981) . As in 
the case of 3C 454.3 discussed bv lVillata etail (120071) . these oc- 
currences might be due to (temporary) different orientations of 
the optical and radio emitting regions, harder events being ob- 
served when the optical emitting region is more aligned with the 
line of sight than the radio one, and hence the corresponding ra- 
diation is more beamed. Alternatively, harder and softer events 
might be intrinsically different. In this case harder events would 
correspond to situations where the perturbation vanishes earlier 
in the jet, sometimes being unable to reach the radio emitting 
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Fig. 3. Comparison among the R-band magni- 
tudes (top panel), corrected for the ELISA con- 
tribution, and the radio flux densities (Jy) at dif- 
ferent frequencies. Cubic spline interpolations 
through the 10-day binned light curves in the 
R and 1 mm bands have been drawn in the 
first and second panel, respectively. Similarly, 
the cubic spline interpolation through the 20- 
day binned light curve at 37 GHz has been 
drawn in the other panels, to highlight the time 
delays between the flux variations at 37 GHz 
and those at lower frequencies. In the fourth 
panel the cubic spline interpolation through the 
15-day binned light curve at 22 GHz is also 
shown (cyan line). Minimum (and maximum) 
flux density levels are indicated with horizontal 
dotted lines. 



regions. Regardless of their origin, hard flares in average should 
be shorter, since they cannot fully develop. 

In this framework we could also explain the sharpness of the 
last radio outburst, and the fast contemporaneous flux decrease 
at all wavelengths. They might be due either to a geometric ef- 
fect, i.e. a jet ben ding in the ra dio emitting region (as proposed 
for 3C 454.3; see IVillata et al.ll2007l) . or to an intrinsic pertur- 
bation power off in the same jet region. We speculate that, had 
this shutdown occurred in an even inner part of the jet, we would 
have not observed the outburst in the radio band, i.e. the event 
would have been even harder. 



4. Cross-correlation analysis 

In the previous section we saw that the brightness increase and 
consequent outburst observed in the optical and near-IR bands 
were soon followed by similar events at mm and cm wave- 
lengths. We can of course estimate time delays by calculating the 
time separation between peaks in the different light curves. But 
this procedure neglects the light curves behaviour as a whole, 
and strongly depends on the time sampling. For instance, when 
comparing the /?-band to the 1 mm light curve, we find a lag of 
about 1 month if we consider the December 2006 optical flare, 
while the lag is reduced to a couple of weeks if we consider 
the outburst peak. Hence, in order to obtain a better estimate 



of the average time delays between the flux changes observed 
at different frequencies, we apply the discrete correlation func- 
tion rDCF: lEdelson & Krolikl 19881 iHufnagel & Bregmanl 19921) 
analysis. Since we are analysing a limited time period with long 
data gaps, we use the cubic spline interpolations through the 
binned light curves that had already been shown in Fig. [3] In 
this way we certainly oversmooth the variations hidden by the 
solar conjunction period, but better define the major signals we 
are looking for. 

Cross-correlation between the .R-band and 1 mm, 37 GHz, 
and 22 GHz splines yields the DCF curves displayed in Fig. [5] In 
all cases the distribution of points shows a significant maximum 
(~ 1) at positive time lags r, which indicates a strong correlation 
with variations at the lower frequency lagging behind the higher- 
frequency ones. The delays corresponding to the DCF peaks are 
20, 60, and 70 days for the 1 mm, 37 GHz, and 22 GHz curves, 
respectively. We obtain a better es timate of the time l ag by cal- 
culating the centroid of the DCF (Peterson et all [19981) . which 
is defined as t c = £ r, DCF, / £ DCF,- for all the points i with 
DCF,- > fcDCFpeak. By considering 0.7 < k < 0.8, we obtain 
lags between 22 and 23 days for the R-mm, 56-58 days for the 
R-37 GHz, and 70-72 days for the R-22 GHz cross-correlations. 

However, we have already discussed in the previous section 
that the 22 GHz outburst seems to be affected by a sort of shut- 
down causing the sudden flux decrease. In the absence of this 
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Fig. 4. Comparison among the i^-band magni- 
tudes (top panel), corrected for the ELISA con- 
tribution, and the radio flux densities (Jy) at 37 
and 8 GHz over the last ~ 12.5 years. 



effect, the outburst peak and following decreasing phase would 
likely have been delayed. Hence, the 70-72 day delay indicated 
by the DCF appears as a "compromise" between the more de- 
layed rising phase and the less delayed decaying phase. Indeed, 
by looking at the fourth panel of Fig. [3] we see that the flux den- 
sity increase at 22 GHz followed that at 37 GHz after about 25 
days, the peak was reached 1 week later, and no lag is visible 
during the flux fading. 

In the 14.5, 8, and 5 GHz light curves, the time of the out- 
burst peak cannot be defined because of the inferior sampling, 
but we nonetheless infer that the ascending part of the outburst 
was delayed by ~ 50, ~ 65, and ^ 90 days with respect to the 
rising part of the 37 GHz outburst, while, again, no significant 
delay is visible in the dimming phase. 

The application of the DCF method to auto-correlate the 
composite /?-band historical light curv^l results in a reinforce- 
ment o f the ~ 8.5 year signal already found by Raiter i et al.l 
(I2006al) . 



5. Spectral energy distributions 



One of the major issues raised by the w orks of lJunkkarinen et all 
(2004) and Raiteri et al. I (120051 l2006allbl) is the existence of an 
extra emission component in the SED of AO 0235+164, in addi- 
tion to the synchrotron and inverse-Compton components. This 
component is represented by a bump in the UV and soft X-ray 
energy range, and its nature is not yet clear. In the case of quasar- 
type blazars, a similar feature is clearly visible only during faint 
states, and it is generally ascribed to thermal emission from an 
accretion disc. However, when these sources undergo a flaring 
state, the bump tends to disapp ear, since it is overwh elmed by the 
synchrotron emission (see e.g. iRaiteri et all r2007bl for the case 
of 3C 454.3). In contrast, by comparing SEDs of AO 0235+164 
obtained at various o bserving epochs in w hich X-ray observa- 
tions were performed, Raite ri et al.l (12006a) pointed out that the 
presence of the extra component can also be inferred in bright 



5 This composite light curve has been obtained by considering the 
iJ-band data back to JD = 2447000, and S-band data converted into R- 
band ones (through th e mean colour index B-R = 1.7) before that date 
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Fig. 5. Discrete correlation functions between the cubic spline 
interpolations through the binned 7?-band and 1 mm light curves 
(blue dots), the 7?-band and 37 GHz ones (red diamonds), and 
the 7?-band and 22 GHz ones (green triangles). Vertical strips 
highlight the range of time lags t included between the location 
of the peak and the centroid of the distributions. 



states, and that the component appears to be variable. But they 
could not discriminate between the two most obvious interpreta- 
tions, i.e. thermal emission from an accretion disc and an addi- 
tional synchrotron component from an inner jet region. 

With the aim of shedding some light on this, we se- 
lected epochs in which the source was in different brightness 
states and in which contemporaneous near-IR-to-UV data were 
available. Magnitudes were corrected for both Galactic and 
foreground-galaxy absor ption by follow ing the prescriptio ns of 
lJunkkarinen et al.l (120041) (see Table 5 bv lRiaiteri et al.l2005l) . We 
then conve rted de-reddened m agnitudes into fluxes using the cal- 
ibration bv lBessell et al.l d!998l) in the optical, and the zero-mag 
fluxes for the UVOT ultraviolet bands mentioned in Sect. 2.3. 

We show three of these SEDs in Fig. [(J compared with those 
corresponding to the XMM- Newton pointings o f January and 
August 2004, published in IRaiteri et all (|2006a). We built the 
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Fig. 6. Near-IR-to-UV spectral energy distributions at vari- 
ous epochs: the SED corresponding to the outburst maximum 
(February 19, 2007) is compared to the faint-state SEDs obtained 
during the January and August 2004 XMM-Newton pointings 
dRaiteri et al.L l2006al) and to two intermediate-state SEDs pre- 
ceding the outburst peak. 



brightest SED with data taken on February 19, 2007, at the 
culmination of the outburst; since there were no simultaneous 
JHK data, we inferred near-IR magnitudes by considering the 
R - JHK colour indices corresponding to the brightest near-IR 
level in the light curves (see Fig.[TJ. In the other two cases (SEDs 
of January 31 and February 5, 2007) all data are contemporane- 
ous (taken within few hours). They show intermediate flux lev- 
els, with a noticeable spectral-slope variation: the near-IR fluxes 
are nearly the same, but toward the UV the two SEDs diverge. 

We found that the least-order polynomial that yields a good 
fit to the near-IR-to-optical part of the SEDs is a cubic; we as- 
sume that these fits, which are displayed in Fig. [6] can reasonably 
represent the contribution of the synchrotron component. The 
UV excess is clearly recognizable in all SEDs, and represents 
the extra emission component. We estimated the contribution of 
this component in the UV bands as the difference between the 
total and the extrapolated synchrotron flux densities. 

Figure [7] shows the relationship between the extra- 
component and the synchrotron-component flux densities for the 
three UV bands: a fair linear correlation can be seen (correla- 
tion coefficient r = 0.996,0.995,0.994 for Wl, Ml, W2, re- 
spectively). Notwithstanding the uncertainties affecting our fit- 
ting procedure, we think that this result indicates that this extra 
component can hardly be explained in terms of thermal emis- 
sion from the disc, which is not expected to depend on the syn- 
chrotron radiation, but more likely it is produced by the same 
mechanism that is responsible for the synchrotron emission. 

6. Discussion and conclusions 

The claim by Rait eri et~aT1 d2001) of a possible quasi-periodic 
occurrence of the major radio and optical outbursts of AO 
0235+164 every 5.7 + 0.5 years led to a mobilization of the inter- 
national blazar community to observe the next event, predicted 
to peak in early 2004. A WEBT c ampaign was organ i sed, wi th 
a huge international participation (Raiteri et al., 2005, 2006a b, 
2007a), but the source remained in a faint state and time series 
analysis on the updated light curves suggested a possible longer 
period of about 8-8.5 years, delaying the occurrence of the next 
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Fig. 7. The contribution to the flux density from the extra com- 
ponent versus that from the synchrotron component in the three 
UVOT bands (see text for details). Lines represent linear fits to 
the data. 



outburst toward the 2006-2007 observing se ason. An increase d 
radio activity was registered in 2005-2006 dBach et all 120071) . 
and a major optical outburst was finally observed in late 2006 
- early 2007, about 8.5 years after the previo us major outburst 
peak, thus confirming the lRaiteri et al.l d2006al) prediction. 

We observed the same event simultaneously in the near-IR 
and UV frequency ranges, and then at millimetric and centimet- 
ric wavelengths, with a progressive time delay toward lower fre- 
quencies. We estimated that the 1 mm and 37 GHz outbursts 
lagged behind the /?-band one by about 3 weeks and 2 months, 
respectively. This l atter lag appears a factor ~ 2 longer than pre- 
viously estimated (Raite riet alll2005ll2006al) when considering 
the historical light curves until spring 2005. This may reflect 
some real change either in the jet structure at still unresolved 
scales or in its energetics with respect to the past. Indeed, there 
are sources that exhibit a characteristic beha v ior for many years 
and t hen suddenly change (see e.g. ISmithl 119961 ; lAller et aTl 
1996). Another possibility is that lack of data for long time inter- 
vals, chiefly in the optical bands because of solar conjunctions, 
affects the results of the DCF analysis. Even in the case of the 
1997-1998 outburst, which was i ntensively observed thanks to 
the efforts of the just-born WEBT dRaiteri et all 1200 ll) . the cen- 
tral part of the event was missed in the optical, and the dimming 
phase was poorly sampled at 37 GHz. As for the 2006-2007 out- 
burst, notwithstanding the exceptional sampling, we cannot rule 
out that solar conjunction hid a further optical peak. This would 
shorten the radio lags. 

In any case, the delays estimated in the present paper ap- 
pear rather short if compared, for exa mple, with the app roxi- 
mately twice-longer lags estimated by IVillata et al.l d2007l) for 
the quasar-type blazar 3C 454.3 in cor respondence of its 2 004- 
2006 exce p tional outburst, and by IVillata et all d2004l) and 
iBach et aTl (120061) for BL Lacertae. The short tim e delays in 
AO 0235 + 1 64 had alread y been noticed in the past dWebb et all 
2000; iRaiteri et all 120011) . and this was one of the reasons why 
microlensing was proposed as a possible e xplanation of the ma- 
jor variability events in this source (see e.g. [Raiteri et all l2007a, 
and references therein). In contrast, when interpreting the multi- 
frequency variabilit y of AO 0235+164 in th e framework of the 
helical jet model by Villata & Raiteril dl999l) . the short time de- 
lays imply that the corresponding jet-emitting regions are only 
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slightly misaligned dOstorero et al.L 120041) . If, besides the geo- 
metrical effect, the outburst is also produced by energetic pro- 
cesses inside the jet, short time lags imply that the emitting re- 
gions are also closeby. 

At lower (< 37 GHz) radio frequencies, the outburst be- 
haviour changes: while the rising phase is progressively delayed, 
as expected, the decaying phase is observed simultaneously at 
all wavelengths. In particular, the 22 GHz flux density reaches 
a maximum value similar to the 37 GHz flux density almost at 
the same time and then suddenly decreases. Although the ob- 
served general multifrequency properties of the outburst can be 
explained by models dealing with shocks propagating along an 
inhom ogeneous jet (see e.g. Hug hes et al.Lll989tlValtaoia et al.L 
1992), the emission behaviour during the dimming phase might 
suggest that the scenario is more complicated than what is de- 
picted in these models. Indeed, it seems as if a kind of shutdown 
occurred in the jet between the mm emission region and the 37- 
22 GHz region. We can speculate that either a jet bending or an 
intrinsic power off of the perturbation may produce this effect, 
which could also account for the different flavours of events ob- 
served in the historical light curves of the source: the harder out- 
bursts, which are stronger at the higher frequencies, suggest that 
a shutdown (of any origin) has occurred between the higher- and 
the lower- frequency emitting regions. Understanding whether 
the geometric or intrinsic scenario would be more plausible is 
not easy. One can envisage that if the flux fading is due to a 
jet bending, the perturbation could continue to travel along the 
jet, following a misaligned path, until it could eventually be ob- 
served at radio wavelengths as soon as the jet (helical) path turns 
again toward the line of sight. For those sources, for which the 
Very Long Baseline Interferometry (VLBI) resolution can sepa- 
rate these regions, we should be able to observe a brightening of 
a VLBI knot. 

We notice that our analysis of the multifrequency light 
curves in terms of harder/softer events p resents some s i milari - 
ties with the generalized shock model by Valta ola~et al.l ([l992), 
who distinguish be tween high- and low-peakin g flares. We fi- 
nally mention that [Hagen-Thorn et ail (l2007al) interpreted the 
photometric and polarimetric variations of AO 0235+164 during 
the December 2006 flare as due to the propagation of a transverse 
shock, accompanied by a small change in the viewing angle of 
the jet. 

Another major issue we investigated is the existence of an 
extra emission component mostly contributing in the UV and 
soft X -ray energy range of the source SED. Swift-UVOT ob- 
servations during the culminating phase of the optical outburst 
show that the behaviour of the UV emission strictly follows that 
of the optical one. Moreover, when constructing SEDs with si- 
multaneous near-IR-to-UV data, and separating the synchrotron 
and extra-component contributions, we found that they are corre- 
lated. Although this result is affected by uncertainties due to the 
decomposition procedure, it nevertheless suggests that an inter- 
pretation of the extra component in t erms of radiation fro m the 
accretion disc is rather unlikely (see iRaiteri et all l2006allbl for 
the various possible interpretations). Also the hypothesis of two 
independent synchrotron components appears now inadequate. 
A still viable explanation could be that of anomalous absorption 
of the optical to near-UV emission in a spectrum that would oth- 
erwise be power-law-like from the near-IR to UV. Indeed, exam- 
ples of noticeable intrinsic absorption are found when analysing 
the spectra of broad absorption line (BAL) quasars, where sev- 
eral lines between ~ 1000 and 1500 A (rest frame) can heavily 
reduce the flux in th e corresponding spectral regions (see e.g. 
iTurnshek et al.LlT988b . A further possibility is that the extra com- 



ponent is the result of inverse-Compton scattering of radio pho- 
tons off the relativistic electrons producing the IR-optical emis- 
sion. Indeed, as stated above, the close correlation and short time 
delays between the optical and radio variations suggest that a lot 
of radio photons are available in the optical emitting region. 
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